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ABSTRACT 


A new method of analyzing Met Rocket wind data is described. Modern tidal 
theory and specialized analytical techniques are used to resolve specific 
tidal modes and prevailing components in observed wind data. The objecti e is 
to formulate a representation of the wind v?hich is continuous in both space 
and time. Such a representation allows direct comparison v;ith theory, allows 
the derivation and other quantities such as temperature and pressure which in 
turn may be compared with observed values, and allows the formation of a wind 
model which extends over a broader range of space and time. 

The current results for a single series of measurements indicates that 
significant diurnal tidal modes with wavelengths of 10 and 7 km are present, 
in the date and are resolved by the analytical technique. The application to a 
wider data sample is required before the full potential of method can be 
determined . 
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SECTION I 
INTRODUCTION 

The primary purpose of this contract is the development of a new method 
for the analysis of upper atmospheric wind and temperature data, and the 
application of these methods to data obtained by the standard Meteorological 
Rodket systems. The method embraces the premise that the most useful fora of 
the data is one which allows a continuous representation in both space and 
time. In order to achieve this representation, the wind data must be resolved 
into distinct components, each component having a distinct set of temporal 
and spatial characteristics. Such a resolution has predictive value, in the 
sense that it can predict the winds at a different place and time, and also 
in the sense that it can infer values of other observables, such as temperature 
and pressure. Furthermore, such a resolution allows direct comparison with 
theoretical predictions. 

The analysis utilizes a data sample consisting of a series of wind 
measurements taken over a time interval of 24 (or more) hours, and covering 
a substantial portion of the altitude range from 20 to 80 kilometers. The 
tools of the analysis are (1) the theory of atmospheric motions in the strat- 
osphere and mesosphere, and (2) an analytical procedure specifically suited 
to the Met Rocket data from Wallops Island. The method requires the resolu- 
tion of the observed winds into distinct tidal modes as well as prevailing 
components. Each tidal mode has a characteristic period as well as spatial 
propagation characteristics. The nraber of measurements is always limited. : 
Accordingly, harmonic analysis can resolve the observed winds according to 
their temporal characteristics only imperfectly. Further, there are several 
inodes present v;ith the same temporal period. Thus, a method to include the 
vertical propagation characteristics of each mode was incorporated into the 
analysis. The method removes modes with short wavelengths, and the remaining . 
long-wavelength modes are then defined by use of harmonic analysis. The 
method relies heavily on the general results of tidal theory, and the analytical 
procedure is based on temporal and vertical "filtering" operations in order to 
resolve the observed measurements according to their temporal and spatial char- 
acteristics. 

A brief survey of past observations and analysis of similar data are con- 
tained in Section II of this report. The analytical method is described in 
Section III and the current results are presented in Section IV. 


SECTION II 


PREVIOUSLY REPORTED TIDAL ANALYSIS 

The presence of tidal forces in the atmosphere has long been recognized 
from the observations of surface pressure changes from which it was determined 
also that the driving, force is thermal since the solar tide is much larger 
than the lunar. The semidiurnal component of the observed solar tide has a 
larger amplitude than the diurnal. For many years, this was the principal 
support for the "resonance theory" which assumed that the atmosphere had a 
natural resonance period close to 12 hours. During the past few years, how- 
everi, it has been determined that the motions are actually forced oscillations 
in vfhich the amplitudes are affected by varying efficiencies of the coupling 
to the energy source. It is now generally accepted that the thermally forced 
tides result from absorption of ultraviolet solar insolation by ozone in the 
stratosphere and mesosphere and from infrared absorption by water vapor in 
the troposphere. Contemporary tidal theory and limited observational data, 
indicate that the amplitudes of tidal components increase from a few centi- 
meters per second at the surface to several meters per second around 50 km 
and become the dominant component of the wind field above 100 km with ampli- 
tudes of over 100 meters per second. The development of the theory was 
hindered by the fact of the equations are nonlinear which prohibit a complete 
global solution and also by the lack of observational data to furnish guidance 
for determination of useful approximations. 

The data are difficult to analyze because the atmospheric wind system has 
several components that vary greatly in character and result from different 
sources. These differences appear in both size and time scales. Some have 
actually been observed and identified; others have been proposed as an ex- 
planation of observational results or on purely theoretical grounds. The 
shortest period, smallest size variations are in the latter group. The 
smallest sizes are usually categorized as turbulence and are generally 
attributed to nonlinearity of flow in high shear regions. The next larger 
scale size is often attributed to gravity waves that may result from a variety 
of sources including many kinds of disturbances in the atmosphere. Although 
these two smallest size components have not been unambiguously identified, 
small variations have been reported by many observers and most treatises on 
atmospheric dynamics include them in some form. 

The next larger scale of variations is thermally driven tides with periods 
of 24 hours and its harmonics. Some of these motions have been identified 
and related through the appropriate theory to the energy source that produces 
them; i.e. , absorption of solar insolation by ozone and water vapor. Diurnal 
and semi-diurnal components have been reported in the Met Rocket data and 
several other harmonics have been reported from other sources. 

Planetary scale waves with periods of several days have been observed also 
These are usually attributed to large disturbances in the atmosphere such as 
magnetic storms . 
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3?he largest period variations are seasonal, semiannual, annual, and quasi- 
biennial oscillation. The source of the 2 -year oscillation is not yet known 
but the major seasonal changes are well described by the thermal X'/ind equa- 
tion in which the zonal wind component is defined by the meridional temper- 
ature gradient. • 

The observational data have been subjected to harmonic analyses in order 
to detect and evaluate components with periods of 24 and 12 hours. The data 
have been analyzed from hourly averages over varying periods of several days 
to months and from series of sequential measurements over a day or two. For 
short periods of observation, the results will be affected by the short 
period irregular components such as interval gravity waves or short period 
planetary xcaves. For larger observing periods , planetary waves and seasonal . 
variations produce uncertainties. In all cases, the best fit parameters were 
determined by a purely mathematical process and then compared to theory. 

Early reports by Smith (1960) who found large diurnal and semidiurnal com- 
ponents from 23 observations in Nevada and Johnston Island, and by Lenhard 
(1963) who reported similar components with amplitudes of 10 m sec"l at 60 km 
at Eglin, Florida were compared to the resonance theory of atmospheric tides. 

The observed phases xjere' not compatible with the predictions of that theory. 

Since the introduction of the forced oscillation theory and the acquisi- 
tion of more data, a large number of tidal analyses of Met Rocket data have 
been reported. 

Miers (1965) reported results from three series of rockets launched 
every 2 hours throughout the day from I'Jhite Sands, New Mexico and from Eglin, 
Florida. Diurnal components with amplitudes from 10 to 20 m sec~l and semi- 
diurnal component with amplitudes of 3 to 5 m sec"l were determined. Reed 
et al. (1966) reported diurnal amplitudes of 5 to 10 m sec”l near tha strato- 
pause from hourly values averaged over several months. 

Diurnal oscillations with amplitude of 12 m sec"^ around 56 km x^ere deter- 
mined from 16 observations xvithin 51 hours at l^ite Sands by Beyers and Miers 
(1966). An associated temperature variation of S^C was also reported but 
attempts to relate the wind and temperature variations through the generalized 
thermal x^ind equation xi?ere unsuccessful. 

The presence of a semidiuranl tide with amplitude of 6 m sec“l at 60 km 
was confirmed by Reed (1969), and Beyers and Miers (1968) confirm tidal com- 
ponents in the equatorial stratosphefe from a series of 24 measurements within 
a 2-day period over Ascension Island. 

In summary, it appears that tidal components with periods of 12 and 24 
hours, are consistently present in the middle atmosphere but the relation of 
these components to the variation of temperature and pressure and to the 
appropriate thermal driving forces remains to be determined. The nexi? analytical 
procedure x^;hich is described in Section III is specifically designed to deter- 
mine these relations. 
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SECTION III 


ANALYTICAL METOOD 


A. INTRODUCTION 

The priina.iry puirpose of the analysis of wind data is to obtain a trepire” 
sentation of the winds that is continuous in space and time. In order to 
achieve a useful representation, the wind data must be resolved into distinct 
components, each component having a distinct set of temporal and spatial 
characteristics. Such a resolution has predictive. value, in the sense that it 
can predict the winds at a different place and time, and also in the sense 
that it can infer values of other observables, such as temperature and pres- 
sure. Furthermore, such a resolution allows direct comparison v;ith theoretical 
predictions. 

The analysis suggested in this proposal assumes a data sample consisting 
of a series of v;ind measurements taken over a time interval of 24 (or more) 
hours, and covering a substantial portion of the altitude range from 20 to 80 
kilometers. The tools of the analysis are (1) the theory of atmospheric 
motions in the stratosphere and mesosphere, and (2) an analytical procedure 
specifically suited to the llet Rocket data from Wallops Island. A brief 
account of pertinent theoretical results and of the analytical procedure is 
given belov;. 

B. THEORETICAL RESULTS 

In the altitude region under consideration nonlinear effects can be 
neglected. Accordingly the total wind is represented by a linear super- , 
position of distinct contributions. These contributions are broken down into 
three general categories; (1) large-scale, long-period motions, (2) tidal 
oscillations, and (3) small-scale, short-period motions. The characteristics 
of these three types of motions" are as follows: 

1 . Large-Scale, Long-Period Notions 

Large-scale, long-period motions var’y slowly or not at all during 
the period of observation. In other v;ords , they are nonoscillatory with a 
characteristic time in excess of 24 hours. Their variation v?ith altitude has 
a characteristic length in excess of 10 kilometers. The thermal wind and 
planetary v?aves are included in this category. Although these have been 
studied theoretically, it is not considered feasible at this time to determine 
their properties in detail from data confined to a single observing site. 

2. Tidal Oscillations 

The dominant tidal motions in this region of the atmosphere are 
the diurnal oscillations (v;ith a period of 24 hours) , and the semidiurnal 
oscillations (with a period of 12 hours). No evidence has been presented to 
date of the existence, to any significant extent, of tidal oscillations vrith 
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shorter periods, Perti^icnt features of five diurnal modes and three semidiurnal 
modes are listed in Table 1, which is based on the results of Lindzen (1967, 
1968). All diui-nal modes have the same temporal period (24 hours), but dis- 
trict variations with respect to altitude and latitude. Two types of diurnal 
modes arc listed in Table 1; trapped and propagating. Mathematically, trapped 
and propagating modes are associated with negative and positive, respectively, 
"equivalent depths." ("Equivalent depths" are the eigenvalues of the tidal 
equation.) Physically, for trapped modes, the response of the atmosphere is 
confined to the immediate vicinity of the regions of excitation. Propagating 
modes are so-named because they propagate energy away from the- region excita- 
tion. The term "phase" in Table 1 denotes the hour of maximum. The ttapped 
diurnal modes, whose phase does not vary with altitude, represent motions vjhich. 
attain their maximum value simultaneously at all altitudes. There is evidence 
(Lindzen and Blabe, 1970) that the trapped diurnal modes are not significant 
beloxi7 the mesopause. 

The following characteristics are true of each mode listed in 
Table 1. The amplitude contains a factor which grows with altitude as (p^/p) 
where p is the local density and p^ is the ground density. This is essentially 
an exponential growth. The relative amplitudes of the eastward and northv;ard 
wind components are a function of latitude only, and are independent of alti- 
tude. Finally, the northward component leads the e as tv?ard component by a 
quarter of a period (i.e., the northward component achieves its maximum a 
quarter of a period before the eastward component). All of the above charac- 
teristics are independent of the nature of the exciting source. For the known 
sources of excitation (due to absorption by ozone and water vapor), the relative 
magnitudes of the different modes are known. Mathematical formulations of 
the stated characteristics are presented in the following subsection. 

3. Small-bcale, Short-Period Motions 

Small-scale, short-period oscillations may be representative of 
internal gravity waves or random motions. Since the continuous spectrum of 
allov/ed gravity waves is broad in terras of both periods and wavelengths, there 
are no general characteristics of such motions that are pertinent in the present 
context. • 

C. ANALYT-ICilL PROCEDURE 
1. Previous Analyses 

Previous analyses of Met Rocket data have treated data separately 
at each altitude; further, each wind component was analyzed separately. Basi- 
cally, the method consisted of removing, in some fashion, a "prevailing" wind 
and a "trend", and then dealing with the residual wind as representing the 
diurnal oscillations (or the diurnal and semidiurnal oscillations together). 
Mathematically this procedure expresses the eastward component of the wind as 

u(z,t) - A^(z) cos Wj^t + B^(z) sinw^t 
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Table 1. Characteristics of some tidal inodes 


Period 

(hours) 


Name and Characteristics 


24 

24 

24 

24 

24 

12 

12 

12 


First Trax>ped Diurnal, Phase does not vary with altitude. 
Amplitude varies smoothly with altitude. Probably not sig- 
nificant. 

Second Trapped Diurnal. Phase does not vary with altitude. 
Amplitude varies smoothly with altitude. Probably not sig- 
nificant. 

First Propagating Diurnal. Approximate viavelength: 28 km. 

Dominant diurnal mode. 

Second Propagating Diurnal. Approximate viavelength: 11 km. 

Third Propagating Diurnal. Approximate wavelength: 7 km. 

First Semidiurnal. Approximate v;avelength: 150 km. 

Second Semidiurnal. Approximate wavelength: 54 km. 

Third Semidiurnal. Approximate vjavelength: 34 km. 
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\ + cos + ’& 2 {z) cos w^t (1) 


r + P^(z) + P^(z) t, 

V7here 

u s= eastward component of the wind 
03^ = frequency of diurnal oscillation 
£(?2 = frequency of semidiurnal oscillation 
z = altitude 
t = time 

Pq = prevailing wind 

p^t = "Trend" ' 

and the amplitudes (C) and phases (d) are determined from the equations 

» ^i =^ arc tan (B./A.) i = 1, 2. (2) 


A similar expression pertains to the northward component of the wind. The 
A*s, B's and p's are treated as parameters to be determined by the method of 
least squares from N observations 

• t . , U Ctjq) • / 

Despite the simplicity of this approach^ results fairly consistent with theory 
were obtained by Reed (1966, 1967), Miers (1965), and Beyers and Miers (1966, 
1968). These results are shown in Figures 1 and 2 (after Lindzen, 1967). 

2. Data Sample and Errors 

For an adequate data sample (say nine or mc'e observations over a 
24-hour period), the analysis can be improved by taking into account general 
theoretical results, as vjell as temporal and spatial "filtering" techniques. 
Some of the details of those techniques depend on the temporal .resolution of 
the data, and on the ability to ascertain the magnitude of the errors in the 
data. In general, it is best to acquire the data at evenly spaced intervals 
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Altitujc-(kin) 



Figure 1. Altitude distribution of the amplitude of the solar 
diurnal cotiponent of u at 30° latitude for winter, 
equinoctial and suinirier conditions. Also shown are 
some distributions based on observations. 




Altitude (km) 



Figure 2. Altitude distribution of the phase of the solar 

ditirnal component of u at 30° latitude for v?inter, 
equinoctial and summer conditions. Also shown are 
some distributions based on observations. 



(e.g., at 3^hour intervals for nine observations). A theoretical estimate 
of the errors in the data is usually obtained from the characteristics of the 
measuring apparatus as well as from the procedure utilized to translate 
direct observations (e.g., radar signals) to v?ind values. Ancrther estimate 
can be obtained by averaging the v;ind values over a specified height interval 
(e.g., half a kilometer) and treating the root-mean-square deviation from the 
average as the error figure. Caving to the difficulty of obtaining reliable 
theoretical estimates of errors, the second approach appears preferable. The 
error figure at each altitude determines the vjeight to be attached to the 
wind values at that altitude. Tor simplicity, consideration of error figures 
and weights are omitted from the following discussion. The methods of error 
assessment and data smoothing w’ith the computer codes are contained in 
Appendix B of this report. 

j/ ' 

3. Analytic Representation of the Measurements 

To formulate the technique of temporal and vertical filtering, 
it is assumed that the measured winds can be represented as a linear combina- 
tion of tidal oscillations j a. slowly-varying contribution, gravity waves, and 
noise. Analytically, - 

2 M(n) - 

u(x,0,t) = e>:p(x/2) S 2 cos w^t + (x) sin w^tj 

n=l m=i “ ■ 

+ P^q(x) + P^^(x)t + 


2 H(n) 

v(x,0, t) = exp(x/2) 2 2 Vnm(0) Bnm(x) cos co t 

n=l m=l 


+ P^g(x) + ]T^(x)t + (4) 

where 

eastward component of the wind 
northward component of the wind 
colatitude (3 = 0° at the north pole) 
local time 



A sin ct> t 
run n 
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r 

X = / dz/H(z) “ "reduced'* altitude 

*^o. 

z “ altitude 

H(z) = scale height of altitude z 


( 0 ^ = frequency of oscillation for nth tidal harmonic, n = 1,2 
U (0) “latitudinal variation of mth mode of the nth harmonic for 


nm 


the eastward wind (knov7n) 


V„„(0) = latitudinal variation of rath mode of the nth harmonic for 


nm 


the northward wind (known) 


~ J^sptesentation of long-period eastward wind 

2 ' ' 

+ P^^t + P^ 2 ^ “ representation of long-period northvjard wind 

G = collective label for short-period contributions to eastward wind 
u 

G = collective label for short-period contributions to northward wind 
= noise in the eastvjard wind component 
= noise in the northward wind component 
M(i^) = number of modes for nth tidal harmonic 


Noise, in the present context, denotes a short-lived transient 
that affects only a few'of the measurements. It is distinguished from the 
short -period contributions C, that may be present in all measurements. The 
coefficients Anrn Bnnr <^^®termine the amplitude (Gnm) and phase (6nra) of 
each mode as follows ^ ^ 


2 2 
Cnm -A + B 

nm nm 


6nro = — arc tan (Bnm/Anm) 


(5) 


Note that the growth factor exp (x/2) has been explicitly introduced as a 
coefficient of the tidal components in Equations (3) and (4). Note also that 
the general theoretical relations between the eastvjard and northward com- 
ponents pertaining to each mode have been incorporated into Equations (3) and 
(4) as well. 
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A. Tempotcii "Filtering” . 

"Filtering" usually refers to an operation on a continuous signal 
In the present case the term is used to denote a linear combination of K dis- 
tinct measurements conducted at the times t}^^, k = 1, K. Let w(k) be a 

set of coefficients and consider the combination 

K 

S (x) I w(k) U(x,t, ) (6) 

k=l 


where u(x,t) is the eastward wind measured at time t. The sum S v;ill be free 
of diurnal and semidiurnal contributions provided the coefficients w are chosen 
so as to satisfy the equations 

• ■ ■ ■ K 

E w(k) cos CO t, = 0, n = 1,2 (7) 

k=l ^ 

K 

T. w(k) sin co t. = 0, n = 1,2 (8) 

k=l " 


Then 


K 


K 


K 


■" “ ‘’uO 5 , '’ul 2 , + ’’u 2 , 5 , 

k-I ^ k=l k=l 


K 


K 


(9) 


+ 2 w(k) Gu(t ) + 2 W(k) Nu(t. ) 

k=l k-1 ■ 


Equations (7) and (8) together impose four conditions on the w(k) coefficients , 
leaving K-4 degrees of freedom. (It is assumed that K > 9). Accordingly 
(K-4) linearly independent combinations of the form of equation (9) are 
formed. From these combinations, the three coefficients P P , P ^ are 
determined by a method of least squares, assuming the contr&utxons rrora N and 
G cancel out. This assumption is justified provided the following conditions 
are met: 

1. The coefficients w(k) are of comparable 
magnitude, so that noise present in a fev7 measure- 
ments does not lead to a spurious result. 

2. The combination docs not enhance excessively 
any short-period contribution contained in G^. 

Let CiiQ be the frequency of such a contribution. 

Forming a combination of the form of equation 
(9) enhances (or attenuates, as the case may be 
Gj^j by a factor Fq, where 
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( 10 ) 




\\ 

These conditions are met provided there is sufficient redundancy 


in the data, tor K=9 (measurements) there are K-4=5 equations of the form (9) 
from which the three P coefficients can be determined. 


The first step of the temporal filtering procedure thus yields an 
initial estimate of the long-period v?ind component separately at each altitude. 


^,.9(x) 


ul 


u2' 


determined at each altitude x. 


The second step is the imposition of the requirement that the long- 
period wind component vary slovrly (and smoothly) with altitude. This is equiv- 
alent to vertical filtering that establishes the smooth representations 




k - 1, , K 


( 11 ) 


A 


These are approximations to the long-period contributions to the wind at the 
times of the observations. 


The next step of the temporal filtering procedure is to foi’m the 


quantities 






( 12 ) ' 


and the combination 


K 


T^(x) = 5 a (k) Q,.(x,t,.) 


. k=l 


u 


(13) 


moved. 


Qy is thus the eastvjard vjind with the long-period contributlon re- 
If the requirement 


K K 

5 o: (k) cos cD„t, = 'S, a (k) sin o^tj^ 
:=1 k= 


= 0 


(14) 




is imposed on the coefficients a (k) , then the combination (13) is free of 
semidiurnal contributions. Since there are K coefficients a (k) and two con- 
ditions imposed by equation (lA) , (K-2) linearly independent combinations of 
the form of equation (13) can be obtained. If conditions (1) and (2) mentioned 
above are observed, then the can be considered essentially free of noise 
and short-period contributions and thus contain the diurnal modes only. 
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From equation (3) 


M(l) 

T (x,i) = exp(x/2) 1 

m=l 


K 

+ B (x) 2 Ci:(kii) sin CD t 

B m /} ^ ^ 

k=l ■ . 


K 

A„^(x) X a(k,jO cos CD t. 
nm It— X ^ 


i = Ifj ...» K-2 


( 15 ) 


Using equation (5) this can be written as 


exp(-x/2) T^(x,i) = 


M(l) 


m=l 




(16) 


, K-2 


The same steps applied to the northward wind component result in the equation 


exp(-x/2) T^(x,i) = 


M(l) 


m=l 


( 17 ) 


Z = 1* 111, k-2 


where 


Fj_(i)^ 


t, (f) == — 

1' CD, 


’ K Y ( f 

S Cf(k, ^) cos cOj^tjY I + I S S(k,f) cos CDj^tj, j 
.k=l- / \k=l! . / 


arc tan 


' K \ 

£ a(k, i) sin cDj^tj^ 

. k=l / 


£ o;(k,f) cos cL’-,t, 

(k=i ^ V 


(18) 


( 19 ) 
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Since tlie notation is cumbersome, it seems appropriate to review 
the significance of equations (16) "and (17). Given k arbitrary coefficients 
a there are (k-2) linearly independent ways of choosing these coefficients 
while satisfying the two conditions of equation (14). These (k-2) sets of 
coeff ic • ents are labeled 

j 

a(k,i), K = 1, k, . £= I, k-2. 


There are, therefore (k-2) linearly independent combinations T 
and T^ that can be formed: 

T^(x,;e), T^(x,.g), ^ = l, (k-2). 


For each i, T^(x,f.) and T^,(x,£) represent, respectively, the total eastv/ard 
and northward diurnal wind at /the time t^(£) , enhanced by the factor Fj^(£) , 

The unknov7ns in equations (16) and (17) are the amplitudes C^,^(x) and phases 
6 ijjj(x) , for m =1, M(l) , where H(l) is the number of diurnal modes. 

In an entirely analogous manner the semidiurnal winds are isolated 

as well. 

5. Vertical "Filtering" 

It has been mentioned previously that trapped diurnal modes are 
insignificant below the mesopause. Thus only the three propagating modes 
needj be considered. These have wavelengths -of approximately 28, 11, and 7 
kilometers. Os'.’ing to the very slow variation v;ith altitude of the amplitude 
factqrs CjLjjj(x) , the shorter-wavelength modes can be removed by vertical filter- 
ing. To illustrate the filtering procedure, let us consider the function 

■ f(x) = cos (cot - ~ x), : 

and the integral 

/ ' ■ ,x+a/2 ''\- 
I(a,x) - / dy f(y) 


X . 2 xa , , 2x . 

= ^ sin “Y" cos (cot - ~ x). 


15 




II: a « X/2, then i(a,x) «0, Thus averaging the combinations and of 
equation (16) and (17) over 7 km and 11 km successively will remove essentially 
all but the first propagating diurnal mode multiplied by known factors. The 


result can be written in the compact 

form 



Cij^(x) cos t£(i) - 8j^j^(x) 

= Y^(x,i) 

£ - ly .. . , k“2 

(20) 

Cjj^(x) sin "t (£) - 

= Y^(x,.g) 

i = 1, . . . , k-2 

(21) 


where Yy and Yy are known. Consider now the graphical representations of Yy 
and Yy versus X. It follows that . 

1. Inhere Y has a maximum; 

u 

6^j^(x) = t^(I), C^^(x) =Y^(x,f) 

2. where Y has a minimum: 

u 

^11^^^ ~ aT” * “ “^y 

3. where Y has a maximum: 

4. where Y^(x,i) has a minimum: 

^ll(x) ~ t^(l) +^, Cj^j^(x) = -Y^(x,<g). 


From the 2 (k-2) such graphical representations a sufficient number 
of points can be established in the vertical profiles of C2^j_(x) and d^Xx) 
to allow a smooth analytic representation of these functions. 

After the longest-wavelength mode has been isolated in this fashion, 
it can be subtracted from the total diurnal contributions of equations (16) and 
(17) and the other modes can be investigated as well. 
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6. Temperature and Pressure Investigation 


It is a general result of tidal theory that, for each tidal mode 
separately, the temperature and pre5>sure variation can be determined from the 
wind components. Explicitly, the pressure variation is 


AP„^(x,0,t) 

nm 


4ao3 

~ — ~ ® (0) “ exp(x/2) 

Rgo)^ ^nm^ H 


X - A (x) cos cn t + B (x) sin cu t 
nm' ' n nm' ' n 


and the temperature variation is 


ATn^(x,e,t) = (h) ^^^(0) i exp (x/2) 


v[[h (x) A^^^(x) -d)A^Jx)] cos co^t 


- ThCx) B - H B (x)1 sin CO tl 

L nm nm' 'j n y 


where 


AP = pressure variation due to mth mode of nth harmonic 
nm 


At = temperature variation due to mth mode of nth harmonic 
nm 


a = radius of the earth 


= frequency of rotation of the earth 


R = gas constant 


g = gravitat ional constant 


H ( 9 ) = latitudinal variation (known), 
nm ' ' 




The altitude-dependent functions and are the same ones 
involved in (and determined by) the wind components. Primes in equation (23) 
indicate derivatives with respect to the reduced height. In taking deriva- 
tives, it is always the case that small errors get greatly magnified. Thus, 
in practice one would integrate equation (23) over a small altitude interval 
(1/2 to 2 kilometers, depending on the accuracy of the data) and obtain a 
jnean value for the interval. 

The total tidal contribution to* temperature and pressure varia- 
tions is obtained by summing over all modes 


/"'AP (tidal) = S 2 



At (tidal) - 2 2 T„„. 

g • . KllU 

* n m 

It should be noted that the same filtering procedures developed 
for the winds cani be applied to temperature and pressure measurements of com- 
parable accuracy . 


SECTION IV 


CURRENT RESULTS 

The methods described in Section III were applied to a series of measure- 
ments from Wallops Flight Center in June 1974. After application of the error 
assessment and smoothing procedures contained in Appendix B, the vertical pro- 
files of the northv/ard component of the vjinds for that series is shown in 
Figure 3. It v;as found that the smoothing procedure had to be applied twice 
in succession in order to reduce the random variations to levels which allov;ed 
the extraction of small tidal oscillations. It was also found that a broader 
data base was required to define the vertical profiles of prevailing zonal 
wind and ambient temperature which are required in the analysis. The values 
of these quantities which were obtained from the CIRA 1965 model are given in 
Appendix A of this report. . 

Plots such as Figure 3 are important to the analytical method especially 
during the initial phases. The location of maxima and minima on the vertical 
profiles are used to determine the relative phases of the predominant tidal 
modes. The values of phases which are actually utilized in the analysis are 
computed from tidal theory for the identified modes. Some examples of this 
procedure for flight nos. 5, 6, 7, and 8 of the June series are sho^m in 
Figure 4 through 7. The phase of the 7 km (diurnal) component expressed 
as multiples of m is shown over the observed altitude range. The solid curve 
was obtained from tidal theory utilizing the prevailing temperature profile. 
The points are from visual observation of maxima and minima on the vertical 
wind profiles. The points are generally alternatively the eastw'ard maxima (o) 
the northv:ard maxima (o) , the eastw^ard maxima (0) and the northward minima 
(A).. ■ 

Verification of the amplitudes of the various superimposed oscillations in 
the June data w-as obtained by applying a standard harmonic analysis to that 
series of data. The results are discussed in Appendix D of this report. 

The computer codes for determining the tidal modes with the filtering me- 
thod are contained in Appendix C. Some examples of the results obtained from 
the June 1974 data are shown in Figures 8 , 9, and 10. The phase of a 10 km 
diurnal componeiit which could not be determined visually as v:a$ the 7 km 
component was, however, determined in the. filtering procedure and is shovm 
in Figurd 8 for the first eight measurements in the June series. 

The short x>?avelength residual after the removal of the 7 km and 10 km wave 
length diurnal modes is shown in Figure 9. This remainder strongly resembles 
a 5 km wavelength semidiurnal mode. 

Finally, the very long w^avelength contribution to the northward v;ind is 
showm in Figure 10. 

Application of the program to data from other observational series at dif- 
ferent time has been, initiated but could not be completed in time for-, this 
report. ■ ■ 


SMOOTHCD NORTHWARD WIND, m/s 


Figure 3. Smoothed northward components of the winds from June 1974 












PHASE (MULTIPLES OF ff) 

Figure 7. ' Phase of the 7 km diurnal mode from flight 8 in Figure 3 (see text for 


explanation and meaning of symbols) . 






SHORT WAVE LENGTH COMPONENT m/s 


Figure 9. 


Short v;ave.l.ength compoucuts after reraoval of 7 kiii and 
10 km components. ’ . :• 
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SECTION V 
CONCLUSIONS 


Pifficulties associated with the organization of .the original data and 
the restructuring of the smoothing process were much greater and more time 
consuming than were* initially anticipated. Thus the amount of actual pro- 
cessed data vjas less than initially expected. However, the primary objective 
of the contract was accomplished in that the development of the new analytical 
piethod V 7 as completed. It was demonstrated that .the computer programs produce 
results V7ith the vertical filtering procedure , which are verified by the totally 
separate harmonic analysis process. Assessmeiat; of the full potential of the 
method must await the processing of a much loiiger amount of data from differ- 
.ent times and locations. 


I 
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APPENDIX A ' 

PREVAILING WINDS AND AtlBIENT TEMPERATURE PROFILES 

■ (■ 

Detailed profiles of the ambient temperature and of the prevailing zonal 
wind have been established by interpolation from tables found in the CIRA 1965 
handbook. Values appropriate to the dates and latitudes of the measurements 
of interest v;ere obtained from the tables at 5 kilometer intervals. Using inter- 
polation techniques, we calculated the required profiles at 0.25 kra intervals 
and smoothed the profiles using the procedures described in Appendix B. 

Graphical representations of the ambient temperature and prevailing zonal winds 
are presented in Figures A-1 through A-8. 
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APPENDIX B 


REMOVAL AND ASSESSMENT OF RANDOM ERRORS IN THE .WIND AND TEMPERATURE DATA 
A. INTRODUCTION AND SUI^IARY 

A smoothing technique has been applied to the wind and temperature 
data supplied by the Wallops Flight Center in order to remove random errors. 
The technique is based on the assumption that, over an altitude interval of 
2 kilometers, physically meaningful variations V7ith respect to altitude can 
be represented by a three-term (quadratic) polynomial. Wavelike variations 
with V7avelengths as short as 5 kilometers are totally unaffected by this tech- 
nique. 'At each altitude, the technique yields an improved ("smoothed") value, 
.the probable error associated with the smoothed value, and the random error, 
defined as the difference betvreen the smoothed and original values. The ran- 
dom error? have been analyzed statistically. In general, the random errors 
associated V7ith the wind measiirements are less than 0.5 meter/second up to an 
altitude of about 55 kilometers. Above that altitude, the magnitude of the 
error tends to increase to a value betv7een 3 and 5 meters/second at 80 kilo- 
meters. The random errors associated with temperature measurements are gen- 
erally less than 0.3°G below 60 kilometers and increase to about 0.8°C at 
70 kilometers. 

The smoothing technique, the analysis of random errors, and estimates 
of errors obtained by other observers are presented in the following subsec- 
tions. The computer codes utilized in the smoothing procedure and the deter- 
mination of random errors is reproduced in subsection 5, along with the 
necessary documentation. 


B. METHOD OF SMOOTHING 

The smoothing procedure is first outlined in general form, and the 
specific form appropriate to the current data is then presented. 


Let Y]^ = the observed value of a physical quantity Y (e.g. temperature) 
corresponding to the value z^^ of the independent variable 
z (e.g. altitude) 

Sj. = the probable error in Yj^ (B-l) 



M 

I 

m=l 


C 

m 



= a suitable analytic representation of Y, 


where f , = f (z, ) , m = 1, M 

mk- m k ■ : 

= the value of the expansion function f at z, 

m k 


Cjjj = expansion coefficient 

M = munber of terms in the analytic representation. 


Define . (B-2) 

I ' " s' 

: . - i ' ' ' , 

' ■ I ■ ■ 

>k " (sX') ^ 

'where K = number of observations of the quantity Y 


The optimum value of the coefficients Cjp_ is obtained by minimizing Q 
xdLtli respect to the variations of, the Cjjj (the least squares procedure): 




Corresponding to the errors in the measured quantities there 
are probable errors in the coefficients Cj^i. It can be shown that these errors 
are represented by the matrix: 

(H-\j • (B-7) 


where the quantity on the left represents the correlation of the probable 
errors and tCj^ , Further, the most probable value of Q is given by: 

Q* = (K - M) ‘ S^. (B-8) 

Finally the most probable error in is given by: 



M M 

E E AC AC„ 




(B-9) 


M M 

= E E (H 
m=l £=1 






If, as in the present instance, the uncertainties Sj- in the measure- 
ments are unknovm, a somewhat different procedure must be adopted. For 
the data under consideration, we assume that the uncertainties in the measure- 
ments are equal over the interval involved. Thus: 

= S^, = 1, k = 1, K (&-10) 


and using (B— 5) and (B— 8) : 


S^‘= E D^/ (K - M) . (B-11) 

For the present data, the analytic representation consists of a three- 
term polynomial v?ith respect to the altitude . The smoothing procedure involves 
nine data points (or equivalently, 2 kilometers) centered at the data point in 


question. This choice of smoothing interval allov;s most of the noise to be 
removed v7ithout affecting wavelike components v;ith vjavelengths as short as 
5 kilometers, vrnen applied to simulated, error-free data resulting from the 
superimposition, of a 5-kilometer v^ave and an 11-kilometer v;ave, the method 
introduced less than 1 percent error, figure B-1 illustrates the effect of 
smoothing on a data sample with particularly severe random errors. Note the 
negative correlation betvreen the eastv^ard and norUhv/ard components. 

c. aI'Ihlysis of random errors 

The random errors established during the smoothing procedure have been 
analyzed statistically. They are presented graphically in Figures B-2 through 
B-5. Figure B-2 pertains to the random errors in the temperature data, plotted 
separately for each series as v;ell as cumulatively for all the data from Wallops 
Island. In each instance the plotted quantity is a 2-kilometer average of the 
root -mean- square va.lue of the random error D as defined by equation (B-2) . 

0v7ing to the nature of the temperature measurement, the random errors are small, 
as expected. Below 60 kilometers the errors are generally less than 0.3°C, 
although the Rourou data gives errors closer to 0.5°G. Above 60 kilometers the 
error increases to a value of about 1°C at 70 kilometers. 

Figure B-3 displays the random errors for the cumulative v;ind data 
from VJallops Island. Belov? 55 kilometers the errors vary betv?een 0.5 and 
1.0 meter/second for both the eastv?ard and northvjard v?ind components. Above 
55 kilometers, errors increase to a value of about 5 meters /second at 80 kilo- 
meters. The dashed curves represent the uncertainty in the smoothed v;ind values 
These are generally about half as big as the random error at the same altitude. 
The apparently oscillatory nature of the random err'or in the eastv.''ard component 
above 55 kilometers is not due to any physical process. It is, rather, a for- 
tuitous result of the superposition of all results from Wallops Island. This 
can be seen more clearly in Figure B-4, which displays the correlated error, 
as defined as: 

D(E) • D(N) 

where D(E) = Random error in eastward v;ind, 

D(N) = Random error in northv?ard v7ind. 

Figure B-4 reveals that above 55 km the eastv?ard and northvTard errors 
are strongly correlated, possibly ov7ing to the experimental difficulty of 
accurately establishing a direction of observation at high altitudes. 

Figure B-5 is analogous to Figure B-3 but pertains to the data from 
Kourou..-, .. 

In suKunary, it can ba said that random temperature errors are insig- 
nifjxant belo'w about 65 kilonieters, v.-hereas random V7ind errors are small belov; 
about 55 kilometers. It should be noted that the technique utilized in the 
current study does not detect systematic errors. . 
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altitude:, km 


Ca) EASTWARD V/IND (M/S) 


(b) NORTHWARD WIND (M/S) 


(c) CORRELATED ERROR (M/S) 


Random Error in the Wind Data From Wallops Island as a Function of Altitude, (a) For 
the Eastward Component. (b) For the Northward Component. (c) For the Correlated 
Error. Solid Curves Depict Error in Original Data. Dashed Curves Depict Uncertainty 
Associated with Smoothed Data. 
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SQUARE OF CORRELATED WIND * ERROR 


Figure B-4. 


The Square of the Correlated Wind Error as a Function of Altitude, Plotted 
Separately for Each Series (see text for definition) . Curves Have Been 
Drawn Only Where the Error Differs Appreciably from Zero. 
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D. OTHER ERROR ESTItlATES 

I Some of the estimates of error associated with standard Met rocket 
Tneasuiiements of wind and temperature are listed in Table B-1. The general 
eetimates were determined by the authors from the various factors involved 
in the measurements. The observed differences \-7ere obtained from comparison 
of tV 70 or more measurements. Susko and Vaughan, tracked Jimsphere wind 
sensors v/ith tv;o different FPS-16 radars and compared the results. Miller, 
Wolf, and Finger, compared the results of simultaneous temperature measure- 
ments made with a U.S. and a Japanese sensor. Miller and Shmidlin compared 
the data obtained from pairs of identical rocket systems launched 5 minutes 
apart, along the same trajectory. Finally, the uncertainty produced by the 
angular^ error of the tracking radar is estimated for data taken at altitude 
increments of 250 meters. 

E. THE COMPUTER CODES USED FOR SMOOTHS THE DATA 

This subsection contains a listing of the computer codes utilized to 
obtain smoothed values of the wind and temperature data, as well as to calcu- 
late the random errors in the ’original data. The smoothing procedure utilizes 
tv70 subroutines, SMOOTH and MISCELZ, together with a calling sequence. 

MISCELZ is called once to establish parameters common to the smoothing of all 
of the quantities of interest. Subroutine SMOOTH is called successively to 
accomplish various tasks, as indicated by the calling sequence. Each sub- 
routine is listed below, together v?ith explanations identifying program vari- 
ables with the mathematical quantities defined in subsection 2 this appendix. 

F. ■ subroutirt: MISCELZ 

Subroutine MISCELZ is listed in Table B-2. This subroutine defines 
the four sets of parameters required to carry out the smoothing procedure. 

The program parameters and their relation to s 5 rmbols previously defined in 
subsection 2 are as follows; 

HI is the inverse of the matrix H as defined in the 

equation preceding (B-6) 

the analytic (polynomial) expansion functions 
defined in (B-1) 

G(M,N) are the auxiliary functions defined by (B-6) and 
the equation that follows it 

ESQ(K) is t he quantity on the right side of (B-9) without 
the S^. 

G. SUBROUTINE SMOOTH 

Subroutine SMOOTH is listed in Table B-3. This subroutine is designed 
to produce smoothed values of the two v7ind components and the temperature. 

It v;as found desirable to do a double smoothing of the original data. In 
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Table B-2. 



SUbROl'TlilE l\:bC“L2 

c cof;TAirs i-.isc;n.LA;:E3^-3 uata 

ccr- r.Ct;/ra SC2/l;l (5,3) iC^ (3i9) ,FSQ(9 ) 

LATA HI/ ’z.bbAlE-Ol, 0» _ ,-3.A652t-01 v 

0» 1 2*6667t*’0l, C* ' 

2 -3.^632E-0l, 0, » y.3I17£-01/ 


LO ]AC. ;:=1,9 

LO 110 f' = 2,3 
110 L = (0.2Ln(l-.-L) 

uO 130 11=1,3 

cc=o.c 

r L=li3 

120 CC = GC+KI (L:,M)r:F(L,iJ) 
130 &(H,f.)=cc : 

140 COraii^UE 
C-~— - 

LO IbO j=l,9 

cc=o,a 



LG l&A M=l,3' 

LD luA :; = 1,3 

184 CC =CC T H I ( I i , L ) { !■' 1 J ) { iv 1 J ) 

185 ESC(J)=CC 
KtTUKH 
END 
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Table B-3. 


SOERU-TImE SrOcTH ('I f I TLK » W.MAX , F A ♦ FB » ES »U5 , E 1 ) 

tIVEiv FA(V.)» }: = ’lMi,--»:.MAX AT 0,25 i'vi-; INTERVALS 

SE.CaiF'S f‘A URi’.'o TiiKuE-'iEA.'-i FuLY.MOMI AL FIT OV'EU TV.Tj l^lLOfiETEHS 
(9 P0H.T5). FI =b:-:0CTnEJ v'ALUEE OF FA ,, 

FOR rTr_F. = 2 calculates Ri-.GuAbLE FRRUf S L'b , ABSULUTi: ERRORS D5' 

FOR I =3 AI.D I T hR = 2 » 3 C ^l.CUL AT E5 F 1 = T ER I VAT 1 VE 0 F FA 

COKFOR/DATAS/VVT (3^361) ,Sr.(3,36l) ,F.SM(3,361) fDSf'.(3r361) ,0111(361 ) 

1 ,D2R<361) 

CC?^r-(;i /l!lSC2/Rl (3,3) ,F (3,9) ,C(3,9) ,ESO(9) 

OlF EL’SlC;; FA (361) ,FlM36l) ,cs ( 56 1 ) , D5 ( 36 1 ) , Fi ( 36 1 ) , C ( 3 ) 

IF I TIALIZE CUT rUT quantities 
LO 19L| • > — ’ li '• I N , f InrtX 
FB (R) =99V.O 
:ES(R) =999.0 
Db(i:) =999.0 
Fl(l>) =999.0 
190 CORTl'- UE 

- JKA^ = 9 : • ■ . 

jF‘1 r = b 
URF =JAAX-3 
JG0=I ■ 

Ji.Mf.= l 

jap; x=jf ax/2+1 

. ■f;LC=L!-.lM-l - 
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Tabl6 B-3 (continued). 

c b^'>co^.l r(..i-:ib AT a 

2CO i:LC-f.LC;+l 

(itiUf AX.. I 

If- (f.Al ,L1 ..’iiiAX) 00 TO 210 

Jf :!•' I ! = Ji-'A X/ 2 I- 1 

Jf';i'.AX = J:--AX 

JOG = 3 

C Mf'L' rXPA.'iGl6f COEKKICli.fnS 

210 GO 23C-. ?!=1 ,3 
GG=0.0 

GO 220 JP=liJMAX 
0 = f:L0 + jr:-l 

220 ct:=CG+G(r.iJr.)i'>i-A(r;) 

230 G(n=CC 

IF ( n t'R,f.!*v2) 00 TO 251 
c CALCOLATr: GFPCO.-S 

i>UM = 0.|C 

GO ?.5C- Jfi=l,Jf.AX 
f. = M.G+jr<-l ^ 
c.e=o,c 

DO 2A0 i;=l i5 
2AG C.C = CC + C (.;) xF (I-., JM) 

GC = CC-V VT ( I ?i’:) 

2tC bU^■=GL.^ +CC-K<i-2 
251 COFTK'UL 

00 TC (270, 260^270) ,JoC 
C J0C = 2 {r- lOOLE I('.TEHv.''.LS) 

260 K=r:Lc+JFA>:/2 
Foc:) =c ( 1 ) 

FKN)=C(2) • 

IF (ITFF.K'd.Z) GO TO 200 
us(P) (1 )-vvT ( no) 

LS (T') =GC.AT (50f ^ESO ( J i lO) /ON;^) 

GO TO 2C0 

C J0G=l,3 {Fir<ST,LAST INTERVALS) 

270 GO 2Fb JN = Ji!i;ri;» JiO-iAX 
K=NL0+JN-1 
CC=0.0 

DO 280 0=1 »3 A 

280 GC=GCtG (!•',) ivF(F»JN) 

FB(N)=CC 

F I CP) =C (2) +2.0:«C (3) , JN) 

IF ( I TE F- , NE , 2 ) GO T J 235 
DS (i!) =CC-WT { I vK) ' 

ES (R) =50r.T (SuNkESQ ( J. l) /Jl.M) 

285 GOATIAUE 

IF ( JGC,Ea.3) GO TO 290 
jGO-2 
GO TO 200 
290 GUNTILUF 
HETLiKN 
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Other vrords, the smoothed values obtained on the first iteration were used as 
the input data to obtain a nev? set of smoothed values. Errors are defined as 
the differ*ence betv7een the original data and the doubly smoothed data. For 
temperature data, the first and second derivatives v;ith respect to altitude 
V 7 ere evaluated as v7ell. Initial values for the first derivative were obtained 
during the secor.d smoothing of the temperature. These initial values v;ere 
then smoothed to obtain the final set of first derivatives and the initial 
values of the second derivative. Filially the second derivative of the tem- 
perature v?as smoothed. 

The subroutine arguments have the follov7ing significance: 

I = index specifying quantity being smoothed 
= 1 for eastv;ard wind 
= 2 for northward viind 
= 3 for temperature 

ITER = Index specifying iteration number 
= 1 for first smoothing 

= 2,3 for subsequent smoothings and derivative calculations 

= Limits on height indices. These are determined by the 
availability of data subject to the restriction that a 
0.25 kilometer grid is assumed. 

FA = Vector containing values to be smoothed 

PB == Vector containing smoothed values 

DS = Vector containing differences between smoothed values 
and original data (as defined by (B-2)) 

ES = Vector containing probable error in smoothed data (as 
defined by equation (B-9)) 

FI = Vector containing derivative of smoothed quantity. 

Within the main body of the subroutine the inde:-: JGO serves to select 
the procedure to be used at the beginning of the altitude range (JGO = 1), 
the end (JGO = 3) and the middle (JGO = 2). This distinction is necessary 
because smoothed values at each altitude are obtained by considering the 
four altitudes that precede it and’ the four altitudes that follow it, except 
at the beginning and end of the altitude range. 
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1. Calling Sequence . 

The calling sequence used to accomplish the tasUs described previously 
is evident from the listing provided in Table B-4. It is assumed that it is 
preceded by CALL MISCELZ which establishes the necessary parameters, as well 
as by the definition of the index I v/hich serves to define the quantity to 
be smoothed. The variables that appear in the listing have the following 
significance. • 

WT (I,N) is a matrix that contains the eastward and northward 
winds and the temperature. 

Sli (I,K) = Matrix containing final set of smoothed values for 

quantity I. 

ESM (I,N) - Matrix containing probable error in smoothed value 

. DSM (I,N) = Matrix containing random error in original data 

DIM (N) = Vector containing first temperature derivative 

D2M (M) = Vector containing second temperature derivative 


no n n n 


::l 


Table B-4. Calling Sequence 


C DLFIi [■ (.UA?;TITY TC u-L 5. OOThr.t) - 

uo 110 

110' f-b (f) -vvi ( I W.) 

C FIH5T bh00ri'l!T 0[-: eUA.vi ITY I 

1TI.A-] 

CALL rrcoTH (I ,it::h,i;vu •:,i-.i;AA»pr,,PK,LS,D5if-i) 

bLCtul’ bLuCTMOG C;F LUA,.TiHY I (LIT! I i^xOLADLL ERROR E5?AD50LUTE 
LIFFLRQRCE Ds. FIRST DERIVATIVE H FOR 1=3) 

m.r; = 2 • 

CALL SFCUTH { I , I TLR , U ; ! ,r-.('-.AX , FR , FS , LS ^ DS , Fl .) 

SAVE CALCULATEL 0UA.<rniE5 
I Lo 120 L=f:sii:,( ;'Ax 
S(-;d iF) =F5(M 
, E5K{ I ,r;) =E5(i.>; 

-::12C OSIW UF) =:[:5(A) V 

IF { I .LT;,3) goto KU . 

FOR T f'f' P r.ix A 1 1 )t< r (1-3) 0 ;LY : 

SHCGTF FIRST TERRERaIORE DERIVATIVE, FIOD-SECOND DERIVATIVE 
ITER=3 

CALL S.f'.CO I li {T 5 I TF.R > I J 1 iRT' AX » F 1 ^ L 1 1 ES t DS 5 Fi' ) 

C SF.CCTR SECOKD TEFPERATOAE DERIVATIVE 

I 7 £R=A 

CALL SfiOCTH { I aTER^iiHlRfR.'IAXiFRfDZF.tES^DStFi) 

140 CORTIf-.UE 







APPEKT)IX C 


CO>D?UTER CODES FOR DETEI^MISATIONS OF TIDAL MODES 

This program resolves short-wavelength modes from vertical wind profiles 
The following subroutines are required: 

1 . EPARV - Reads parameters for vertical analysis. 

2. RFILEV - Reads data from smoothed vertical profile. 

3. RPRE\^ - Reads prevailing w'ind, temperature and 

temperature derivatives. 

4. ZFUN3 - Defines height as a function of height index. 

I' . ' 

5. VFUKS - Calculates .expansion functions for 

vertical analysis. 

6. VFilter - Extracts short-wavelength modes. 


000003 

OOCC03 

000003 

000003 

0000:03 

000003 

000003 

000003 

OOOOG/. 

000005 

0COC06 

000007 

0000 II 
000012 
000016 
000016 
000030 
OC0030 
00C036 
000036 
0000^0 
000056 
000056 
000065 


000065 


000075 

000107 
0001 10 
000112 
000122 
000123 
000125 
000127 

000130 

000131 

000133 
000I3A 
000135 
0001 AO 


PP.OCKAK. TMODESC lNPUT,0'JTPUT»TAp£l = inr’UT,TAPt 2 = 0 UTPUT«TAPElU 
1 TAPL12) 

C 

C Rf:5^LVf:S SttOPT-WAVClENOTH f.QDES EE-OM VlRTICAL PROFILES 

C RFOU I KES RPARV ,KF 1 LEV , VFUfJS , ZFUN , ARRAY , M i t;V , PHASE 

C 

COr.MOr/iOOf;V/fJIT,NCT,KTAPC,LTAPE,PI ! 

CO^’MC/f./HFALiS/ I SLR ,L , L5YAX , TCL ( 7 ) i FLAT , 1 LCHG, TLOC ( 1 a) , ,NHR ( 1A ) , 

, I L'i-' I ( lA) ,iluAY (16) ( 16) ,lVE ( 16 ) if (01 ( 16 ) , fi TP ( 16) , vehicle (16) 

COMr'ON/OATAS/VVT( 3,260) ,SM(3,260) ,EbH(3,260) , OSH (3, 2 60) ,Dlf’U260) 

1 <U2.'-i (260) .Af.N (3i 16 ) ,M-.X (5 1 16 j 

COMMOK/PKEVS/PKEV (3.260) .’.MAX 

COM!:0(l/PARV/KP0L,Zu,SCALE,f;HR,NH(3) .MHO) ,HHR (3) ,NOPT»hlNT 

CaMFOr:/CoDES/n-PR:U,iFTAPE 

COMI-',0fi/FH.'ls/FKN(7.260) t AL (2 . 260 ) f X I (260) 

C 

MIT = 1 
HOT = 2 
MTAFE=ll 
LTAPE=12 
PI=3, 16159265 

C READ PARA-MeIERS FOR VERTICAL ANALYSIS 
, call RPAf-V , • ■ i 

WRITE (NOT, 100) 

100 FCTMAT (iHi/* PARAMETERS FOR VERTICAL A.NALYSI5#/) 

V.'RITE (NOT, 110)! <P(j!,, ZB, SCALE • ■ 

llO F0".‘'.AT (13, j.; POLYNOMIAL TERMS, ZS = ft ,F5«2 ,* SCALE = A, FA , 2 > 

WRITE CNOT,120) AHR 

120 FORI'.AT (/l3v« HARMONICS«//it K N K HKN*) 

, DO 130 K=l,K.HK . 

130 WRITE (NOT,1AO) K,NH(X),MH(K),HMN(K) 

160 format (3l3,Fa.6) 

V.RITE (NOT, 150) .’•JCPT,KINT ‘ 

150 for, MAT Cr3,« = OPTlO.N, HINT-*,F5,2) 

C HEAD' PRINT /Tape codes 

PIEAD (NiTtlAO) IFPRNT,IFTAPE 

C READ FILE NUMBER LI.MIT5 

READ (NIT, 160) LMlN.LMAXtLS.MAX 

C— — — — 

call rprev 

DO 160" N= 1 ,260 
1613 SM(3,M)=PREV(3,N>-273.0 

L=1 . 

NMN(3,U=r 
flKX (3.L) = 260 
TLOC(L)=0.0 

C , DEFINF expansion functions , 

CALL VFUNS 
DO 600 LL = LMIN,LMA:X 
C READ DATA 

call RFILEV 
call vfilter 
800 continue. 

END 





c 

c 
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000002 

c 

000002 

OOOOIA 

c 
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00002A 
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0C00A3 

000052 

000053 


sunPouTiiji: rparv ; 

p.lad i^akametlihs fop viktical analysis 

2A MAY 75 

COMMO!;/ I QOEV/;j I T ^ .'iOT , MTAPR , LTAPF vP I 
COMKor,/PAPV/l''POL,ZO, SCALE, kHN,KH(3) ,fiH(3) ,HHNO) 
READ rUMBEP OF POLYNOfllAL CCeFF I C I EliPS 
REAP (MIT, 100) kP0L,2B, SCALE 
100 format (I3,2F;3.a) 

READ harmonic FUNCTION DATA 
READ (NIT,120) KHP 
DO 110 Krl,KHR 

llOREAD (NlT,l20) NH(l>>),Ml-i(K),HKN(K) 

120 FORMAT (2l3»F8,A) 

READ OPTION (NCPT--1 FOfN FIXED FITTING INTERVAL' 

RE'D (NIT»100) NOPTiHINT 

RETURN 

END 


ORIGINAL PAGE IS 
OE POOR QUAUIEB 


.NOPTiHINT 


HINT) 


I 


I 


000CC2 

000002 

000002 


000002 

000055 

000055 

000057 

CC0066 

0000C6 

00OC7? 

000075 
0001 1 1 
000125 
OOOlAl 
000143 
000156 
000171 
0C0171 
000173 
00C201 
000201 
000 2 G4 
000212 
000212 
000213 


SUBKCUTlflE f^riUEV . i • 

C READS [TAtA from SrOOThCO VERTICAL FILE ^ 

C 23 FAY 75 

CCF'f'OK/lt OEV/NIT,AOT,MTAPE,LTAPE,f'I 

C0i”.f-0i*> 2*/ I Sr. R 1 L 1 L 5’''A.X 1 TCE { 7 ) L-‘'T » t LCifiCi «.T.LOC ( 1 4 ) 

I f.f-'I (U> ,;;L;AY(14) ,NY0(14) ,KYE(14) ,f;rl (14) ,NT?(14) ,VEHICLEC14) 

COY.f'.Cn/OATAS/VVT (3,26,0) ,5Mt3,2601 , LSi'. I 3 , 260 ) ,D5F(3,260) , DlK (2:60 ) 

1 ,D2.M (26U) ,AY,M3, 14) tl.Y.X (3, 14) 


210 


220 


230 

240 

130 

140 

150 

160 


READ ( i'.T A P E , 2 1 0 ) I 5E R , L , L5MAX , T C 5 , F LAT , FLCMG , TLOC ( L) , NHR ( L ) , 
1 f.Mi (L> ,?;oay(L) ,;p-'0(L) ,nye(U ,rnP(L) ,afi (d .VEHiCLE(L) 
FGra'iAT (3l3,7Al072Fli,6,F9.4,7!3.A9) 

DO 240 1=1,3 

READ (f(TAPE,220).NMIN,NMAX 
FCRFAT (2 i 3) 
liiiRd ,L)=RDI’1 - . 

r.nx ( I ,L) =i;!lAX 
HEAD (.".TAPE, 230) 

(.".tape, 230) 

(.’•‘•TAPE,2 30) 

XE.3) GO TO 
(f(TAP£,2 30) 

(MTAPE,230) 

format dope. 2) 

CC.MTir.-UE 

HEAD (MTAPE, 130) A 
FCEKAT (Ml j 

IF (EOF, MTAPE) 160,140 
V..RITE (MOT, 150) L 

FORMAT EOF NOT THERE AFTER FILE*, 1 3/) 

RETURN' 

END 


READ 
READ 
IF (I 
READ 
READ 


( SM( I vN) ,M=MMIN,NMAX) 
lES.M ( I ,N) ,.\' = ,\M1N,NHAX) 
(VVT d ,N) ,N = NMIN,NiMAX) 
240 

(Dlfi(N) ,N = .NMIN,NM.AX) 
(D2M1N) ,N=N,MI.N,f,.MAX) 


OBIGHAL B 

OF sooR quaifh:. 


000002 

000002 

000002 

000002 
000002 
000002 
C00C03 
00000 ^ 
000012 
000012 
000014 
000017 
00002 1 
000C36 
000035 
000037 
000056 
00CD7A 
000112 
0001 
000115 


I 


1 


j I 


T ' ■ ■ . k 


c 


SUBROUTirJE 

RhADS PRf-VAlLlAG V.-lfJD,TE’-*PEHATUnE ARD TEMPERATURE DERIVATIVES 
CCr.MOfl/IODEV/fJITiMDT ,MTAPE,LTAPEtPl | 

C0?^M0f:/HE/*DS/ISCRtL,LSMAX,TCt><7) , FLAT , FLOMG , TLOC ( 1 A ) ,i'iHR{lA) , 

1 Ml*. I (14> ,MUAYUA) ,i.‘MO(lA} ,;;YLUA) ,f:FI (ia) ,NTP(1A} tVEhlCLEClA) 

CG!'.MOM/DATAS/VVT n,260) ,SM(3,260) ,£SM(3,260) ,D5M (3,260) , DIM (260) 

1 ,D2M(260) ,il?.;M(3ilA) ,.NMX(3»1A) 

COMMOO/PREVS/PREV (3,260) ,NM IN, NMAX 
DIMENSION DUMMY (SO) 

NM'Nr 1 ■ 

fi(-*AX = ?60 . 

READ (NIT, lOO) ISEk 
100 FGRKAI (13) 

DO 150 I = 1,3, 2 
READ (NiTi'lOO) 

IF (I,EQ.3‘) GO TO 120 

read (NIT, llO) (PREV(I,N) ,K = NMIN,NMAX) ,DUHK 
no format (lOIFb.21 
GO TO 1 50 

120 READ (Kl'T,110) DUMMY, (PREV ( I ,N) ,N=NMIN,NMAX) ,BUNK 

READ (N IT, 110) DUMMY, (Dim (N> ,N=NMIN,NMAX) ,BUNK • 

READ (NIT, 110) DUMMY, (D2M{N) ,N=NMIN,NMAX) ,BUNK 
150 continue 
RETURN 
END 



60 


FUfJCTIOM ZFUN(M) 

C DFF‘I(i[:S FitlKjh'T AS A FUNCTION OF HE lOHT- 1 MDpK 

C ZFUf] ( I ) =20,00 KM, ZFUN(2)=20,25 KM, ..., ZrUN (281 ) =90.00 KM 

000003 ZFUN=19.75 + 0,25^«N 

000007 RETUKN 

000007 END 


1 


I 


1 


I 


I 


000002 

000002 

000002 

000002 

00C0G2 

00G002 

00G0O2 

000002 
000006 ■ 
, OOGOll 

000012 

000014 

000017 

000024 

000025 


000045 

000046 

000050 

000051 

000052 

000054 , 

000062 

000064 

0C0P65 

000066 

000070 

000072 

000075 

000076 
000077 
000 LOT 
000106 

000106 

000113 

000115 

000117 

000122 

000124 

00G130 

000136 

000140 

000142 

000144 

000153 


( 


SUr-ROUTINE: VFUN5 

C CALCUlATr (:X»’Ai‘)5I0rN FUi'lCT I (j?* 5 FCf< VfJ<TlCAL ANALYSIS 
'COMf-'Of./lOLVhV/HIT ,NCT ,MTAPLvLlAPE ,FI 

COMMOli /HLAOli/ 1 SLR , 1. , LS.'-'.AX ♦ TCB ( 7 ) i FLAT , FLOf.G, TLOC ( 1 4) ,NHR ( I M , 

1 IK'I (14) vUDAYdO iF-R.'jLW) ,NYL ( 14) -iNF I ( 14 ) ,nTP(14) *VEHIClE(14) 

COr.r.fJ!/DATAS/VVT (3,260) , SR ( 3 , 2 60 ) , ESM ( 3 , 2 60 ) ,D5M ( 3 , 2 60 ) ,DlM(260) 

1 • iD 2 f, (260) ,Y!l'Ui(3,14) ,Nf;x(3*14) 

COnROR/PARV/KPOL , » SCALE ,KHR ,NH ( 3 ) *RH ( 3 ) , URN (3 ) ,N0?T *HINT 
C0l'.?'U)r,/C0Dl':5/IFPHN| , IFTAPE 
COMKCN/FKNS/FKri.(.7,260) ,AL (2^260) 1 XI (260) 

Oir.EFSlCN ALP(3) «0XN(3) ,0XH(3) ,CXR(3) *WAVE(3) 

e 

Kr'.IN-Nl-'.MO^L) _ 

KMAX = IlfiX (2 ,L) 

IF (riPIN.EU.MMAX) GO TO 500 . 

C DEFIfir POLYNOMIAL FUNCTIONS 

DO 100 M=NMIN,NMAX 
FKri(l,M)=l,0 
Z=SCALEft (ZFUN(N)-ZB) 

DO 100 K=2,KP0L ' 

100 FKN(F,N)=Z«A-CK-1) • 

C DEFlNFUAKMUNlt FUNCTIONS 

CHT=2.932E-02 
C6AM=2. 0/7*0 
DZ=0.25 

NMip^roi r . 

C INITIALIZE 

DO no K = 1 ,KHR 
110 AL (K,NHlN) =0.0 
XI (NMIN) =0.0 
LINES=60 
ITER=1 

. 120 N=HNIN-1 
130 N=f! + 1 

IF (N.GT.MMAX) 60 TO 180 
IF (NjEQ.NIllN) GO TO 150 
C FOR Ol.GT.NMIN) REDEFINE QUANTITIES 

■ HH=HNY ■ 

DO no k=i,khr 

140 OXH(K)=OXN(K) • . 

150 continue 

C .CALCULATE NEW SCALE HEIGHT AND ITS DERIVATIVES WRT HEIGHT Z 
HN=CHTf.\5M(3,N) +273,0) 

, - Hl = CHT*^DlM (N) - 

H2=CHT»D2M(N) 

- :::PH = HN+( (HI+CGAM) 

C CALCULATE t‘)EW WAVENUMbERS and WAVELENGTHS 

Do 160 K= 1 ,KHR 

CC=-0.25+PH/HMN(K) 

160 OXN <K ) =S0RT SCO 

IF (ITER. EO, 2) GO TO 200 

C FOR first IITERATION, FIND UNSHIFTED XI, PHASES 

IF(N.FO.NMIN) go to 130 
CC=0.51vD2 

XI (N|=X1 (N-1)+CC»K1.0/HH+1,0/HN) 

DO 170 X=1|RHR 




000155 


170 

AL (K,M) =AL {<,N-1) *CC<K0XH(K) /HH+OXN (K')/HN> 

000177 



GO TO 130 

000177 

C 

180 

IF (ITFH.EO.2) go to 500 

shift XI, PHASES AT END OF FIRST ITERATION 

C00201 



XIM=XI (NMIU) 

00020 A 



DO IF'-. K = 1,KHR 

000205 


18 5 

ALP(K)=AL(K.NMID) 

0002 15 



DO 190 N=NMIN,NMAX ' 

000217 



XI i;\) =Xl (N)-XIM 

000222 



DO 190 K=1,KHR 

000223 


190 

ALCX,I,')=ALP{K).AL(K,N) 

000236 



ITER=2 

000237 

c 


GO TO 120 


c 


SECOND ITERATION ONLY 

000237 


200 

Z=ZFUM(N) 

0002A2 



EX=EXr«0.5«XI(N) ) 

00C2A7 



PHI = 0 , 5* ( RH 1 ♦ H2 2 ) 

00025A 



DO 210 K=-1,KHR' 

000255 



OXR (K) =PH1/ (HMM(K) *OXN (K) <^*2> 

000262 



CC = EX/5GRTrOXM(kJ ) 

000266 



CC=1.0 

000267 



M = KPCL + 2«K-1 • 

000272 



FKN{M,M) =GC^»COS(AL (K,N) ) . 

000302 



M = |-U1 

000303 



FKN(M,n)=CC«SlN(AL(KfN) ) 

000313 



WAVE<K) =2t0^PI*:h;V0XN(K) 

000320 


210 

ALP (K) =AL(K«'N) /PI 

000326 



IF ( IFPRNT,EQ.i ) GO TO 3p0 

000330 



IF (LIf’ES.LT.SO) GO TO 250 

000332 



WRITE (NOT, 2201 I SER , L , TLOC ( L > . ; 

0003A5 


220 

FORMAT (IHI ,^fSERIE5«, 13 (li FL I GHTit , I 3 » TLOC = <t , F7 . 3/ 



1 3X*1HN,5X» 1HZ*5X, 1HH»AX»2HH1 tAXtZHXl »2 (5X»2HOXi5Xt 

2 2HQR,3X,AHWAVE,2X,5HPHASE) ,2 (5X,2HFCt5X,2HFSn 


0003A5 



LINES = 3 

00C3A6 


230 

MfilN'=KPOL + l 

0003-50 



NMAX=XP0L+2«KKR 

000352 



WRITE {NOT,2AO) N , Z ,HN ,H 1 « XK N ) , ( (QXN(K) ,QXR (K) ,WAVE(K) ,ALP{K) ) 
1 K=l,KHR) , (FKN(M,N) ,M=MMIN,MMAX) 

000A22 


2A0 

format (1.HH,I3,AF6,2, (12F7.3) ) 

000 A 22 



LINES = LINES+l- 

000A2A 


300 

CONTINUE ■ * 

O0OA2A 



GO TO 150 

000A25 


500 

RETURN 

OOOA26 



END 
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000013 
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000017 
0CO022 
00002^ 
000037 
OOOOAI 
0000/f5 
000054 
000C62 

000064 

000076 

000076 
C0C104 
0C0i06 
000110 
000111 
000112 
000) 14 

cool 16 
000120 
000122 
000124 
000125 
000126 
000127 

000127 

000133 

000134 

000135 

000137 

000141 

000153 

000155 
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I 




\ I 


I 


C 

c 


c 


c 

c 


SUDROUTinE VPILTtR 

F.XTKAC.TS SllOrU WAVCLCNGtH flODES 
COI-*.MOfl/ !n.)pV/ri I T ,M 01 »K,rAPi: ,LTAPF ^ P I 

COP.KCf i/H KA D-5 / I 5 fP , L , L 5.''iA X , T C13 ( 7 ) , f L A T , P L Of .G , T L OC ( 1 4 ) , HR ( 1 4 ) , 

1 fff’.I (14) ,flDAY(14) ,r;f‘0(L41 ,NVE(1 a) ,KFI (14) ,NTP(14) »VEHICL£(14) 

COIlMOf</DATAS/ VVT ( 3 , 260 ) , B'<-H 3 , 260 ) i EG.M ( 3 , 260 ) , DSM ( 3 , 260 ) , DlM ( 260 ) 

I v02M ( 2 60 ) , f.r.t-; ( 3 , 1 4 ) t ftpx { 3 » 1 4 ) 

CCf-rOfi/PARV/KPOL ♦40 t SCALE! , KriR, N’K ( 3 ) ♦ MH ( 3 ) , HMN { 3 ) » .NQPT * H INT 
COMI'.Ci! /C0i>E5/ I PPRM , I FTAPE 
COr.PCrf/FKNS/i'KP (7»260) i AL (2»260) iXI (260) 

COMMON/PREVS/PREV ( 3 I 260 ) , NM I N , RMAX 

OinPf-SlOM H ( 10, 10) ,KI U0*10) ,U(10) ,C( 10) »SV(100) »AMP(3) »PH(3) 
DIP.PMSIOJ) YC(fn) »ER(81) 

DiPEMSIOiJ y!<260)‘,*>M260) 

DO 4C0 I=l,3 
( I ,L) 

NMAXrf.flXd ,,L) i ! 

IF (fiP.in.E0,H(-1AX) GO TO 400 
IF (I.E0.2) GO TO 110 

CCrO.O " • • 

IF U.EQ.3) CC = 273.0 

PO 100 N = NiMIN,.NPiAX . 

100 s;;(i,r:)=sH(i,N)-pREVn,M)+cc 

110 DO 120 N=nMlil,}ir-iAX 
Y(f;) =r,p, (i,N) 

IF <ESM(I ,N) .LT.1.0) ESM(l,N)=l,0 
W(f;) = 1.0/ (ESt-K I ,N)*ff2) 

120 COMTIKUE 


WRITE U10T,130) ISER,L,I 

130 FORMAT (lHl/3l4,i» = SERIES, FLIGHT, OUANT I TYx-/3X , IHN ,7X , IHZ , 6X , 

1 2hYC.7X,1HY,6X,2hEr,5X,3HSUM.6X,2HYP,2(5X,3HAM,P,6X,2hPH) /) 

JMAX = 4« (HlfiTrO.OD + l . 


KMAX=RP0L*2^K.HR 

DfiP.sJfTAX- ( ai'VAX) 

JG0=1 

JllMIfi=l 

JfiM.AX = JMAX72-*.l 
NLO = Mf'lN-l 
200 ML.0=NL0 + 1 

NHI=f!LC + jMAX-l 

IF (MIII.LT.RMAX) GO TO 210 

JNMIN=JWAX/2+l 

jmmax=jmax 

JG0=3 



210 COIITIMUE 

DeFIKE/REDeFIK’E MATRIX H, VECTOR U 
IF (PLO.GT.MMTN) go to 250 
FIRST DeTIMITIOR OF H,U 


DO 240 K=1,KHAX 
GC=0.0 , 

Do 220 JTi=:i,JMAX 
N = riL0 + jN-l 


2 2 0 c C = C c + w ( n ) « F K n { K , N ) Y ( N ) 
U(K) =CC 

DO 240 K=i,XMAX 


6A 


000156 
000157 
COM 61 
000163 
000177 
ooo?ib 
000210 

000210 
000212 
00021 A 
0002 15 
0C0235 
000236 
000270 

000270 

000272 

000273 

00C303 

000306 

0C03D7 

000317 

000321 

000322 

000324 

000336 

000342 

000343 

000344 

000346 

000347 

000351 

000362 

0X10365 

00C370 

OC0376 

C00AO2 

000411 

00C413 

000415 

000417 

000422 

000423 

000424 

000435 

000437 

000442 

00045 2 

000->63 

000466 

000527 

000527 

000532 

000534 


I 


1 


f 


1 


CCrO.O 

DO 230 JI1=1*JHAX 

M=ru.o+jD-i 

230 CC = CC*WU1)*FKN{K,N)#FKN{MiN) 

240 H (K ,(•'.) =CC 
GO TO 270 
250 COIlTirUE 

C HEirtf IHE H,U • 

NF=NLr-l 
NL = flllI 

DO 260 K = 1,:<MAX 

U ( K ) =U ( K ) *W ( iNL ) tf-FKN (K » NL ) *Y (NL ) -W (Nr ) -(tFKN ( K , NF ) Y CNF) 

DO 260 M=1 ,14MaX 

260 H (K,M) =H (K.,1-1) *w <iiLKFKN (K,NL) *FKN (MtNU-W(NF HtFKN CK ♦NF ) ♦fFKN (M i NF ) 
270 CONTIMUE 
C liMVliftT matrix H 

KK = K‘-iAX>i'RilAX 
MGDE=2 

CALL array (MODE, :<MAX,KMAX,,KKao» 10 »SV,H) 
call niNV(SV,KHAX,RK,D)i 
MOOEi: 1 ' 

CALL ARRAY(MCDE»DMAX,XK-AX,KK,10tlO,SV,hI) * • 

C FIMD FXPAiMSIOM COEFFICIENTS 

DO 290 K=1,KMAX . ’ 

CCrO.O 

DO 2F0 M = 1v<MAX 
280 CC = CC + IJ(M)*HI (M,K) 

290 C(K)=CC 

c calculate errors 

SUM=0.0 . 

DO 310 JN=liJMAX 
K=NL0+JN-1 

CC=0.0 

DO 300 K = l ,Xt'lAX 
300 CC=CC+C(K)*FkN(K,N) 

YC(JN)=CC 

ER(JM)=CC-Y(N) 

310 SUtl = SUM + W(N)frER(Ji'nJk«2 
5UM=SCRT (SUM/DNM) 

GO TO (330,320,330) ,J60 
320 JnMIf; = JMMAX 

330 DO 370 JN=JNMIN,JNMAX • 

M = r!LO,.jN_l 

Z = ZFU(.(M) ■ 

YP-0.0 

DO 340 K=r,KPOL 
340 YP=YP + C(K)itFI<LN(K,N) 

DO 350 K = 1 ,KHR 

KK = KPCL+2<»K-1 ! 

AMP (K) =5Q!n (C (KK) ^;-*2 + C (KK + 1) «*2) 

350 PH(K)=PHASE!C(KR*1) ,C(KK))/PI 
YY=YC(JN)tYP 

WRITE (MOT, 360) N,Z»YC(JN) , Y (N ) , ER ( JN 1 » SUM, YP , ( ( AMP (K ) ,PH (K ) ) ,K = 1 
1 ,KHR) ,YY>2 

360 format (14, 12F8. 2) 

370 CONTINUE 

IF (JGO.EQ.a) GO TO 400 

JG0=2 • 
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0CC535 oO TO 200 

000535 400 COWTIf.UE 

000537 return 

000540 END 








APPEM)IX D 

HAPJIONIC AIJALYSIS OF WIND DATA 

Standard hartronic analysis techniques have been applied to the data from 
the June 1974 series at Wallops Island. The results are of a preliminary 
nature, in that the vertical resolution techniques have not been applied. 
Nevertheless, the preliminary results are sufficient to demonstrate the valid-' 
ity of the elements of tidal theory that have been previously mentioned. Of 
particular interest in the present context is the fact that the eastv;ard (X) 
V7ind component lags the northward (Y) V7ind component by a quarter of a period. 
Let the temporal behavior of the wind components be expressed analytically as: 

X = C coswt + S sincPt (D“l) 

X X 

Y = C cosopt + S siiiOJt. 

■ y y 

Then it follov?s from tidal theory that: 

C = K • S S = -K • C , (D-2) 

y X y X 

where K is a constant determined by the latitude. This result holds for each 
tidal mode individually. 

I 

The correlation implied by Equation (D-2) is evident in Figures D-1 and 
D-2, Figure 1 displays the coefficients Cj-j S^., Gy, and Sy for the diurnal 
v;ind component, as determined by harmonic analysis , It can be seen that the 
correlation between C.^ and -Sy, and between Cy and is good, even though 
the diurnal component is made up of several modes. In Figure D-2(a), the 
correlation in the semidiurnal component is even more remarkable. Some 
degree of correlation is also evident in Figure D-2(b) for the terdiurnal 
component . 

Figure D -3 displays the prevailing winds for the eastward (curve (a)) and 
northv7ard (curve (c)) v7ir(d components. Also plotted at 5 kilometer intervals 
are values of the eastward prevailing wind for July 1 at the latitude of 
Wallops Island , as listed in the 1965 Cos par Reference Atmosphere. The gen- 
eral agreement betv:een these values and the prevailing wind deduced by harmonic 
analysis is good. The significance of the eastward trend (curve (b)), is un- 
'clear at this time. The same remark applies to the nortiward prevailing wind 
and trend (curves (c) and (d), respectively). 

In sunmiary it can be stated that the preliminary results of the harmonic 
analysis are in general agreement with the_ results of tidal theory. 
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Figure D-2. Harmonic Coefficients for the Semidiurnal and 
Terdiurnal Wind Components. 
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